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Abstract 

We study the agglomerative percolation (AP) models on the Bethe lattice and the triangular 
cactus to establish the exact mean-field theory for AP. Using the self-consistent simulation method, 
based on the exact self-consistent equation, we directly measure the order parameter P^ and 
average cluster size S. From the measured P^ and S we obtain the critical exponents j3k and 7^ 
for k = 2 and 3. Here f3^ and 7^ are the critical exponents for P^ and S when the growth of 
clusters spontaneously breaks the symmetry of the /c-partite graph [2]. The obtained values 
are /3 2 = 1.79(3), 72 = 0.88(1), /3 3 = 1.35(5), and 73 = 0.94(2). By comparing these values of 
exponents with those for ordinary percolation ((3^ = 1 and 700 = 1) we also find the inequalities 
between the exponents, as (3^ < f3% < fa an d 700 > 73 > 72- These results quantitatively verify 
the conjecture that the AP model belongs to a new universality class if symmetry is broken 
spontaneously, and the new universality class depends on k [Lau et al, Phys. Rev. E 86, 011118 
(2012)] . 
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I. INTRODUCTION 



Percolation transition describes the emergence of large-scale connectivity [lj]. It has been 
extensively studied in various branches of science due to its wide range of applications 
to many phenomena such as sol-gel transition and polymerization, resistor networks, and 
epidemic spreading The first theoretical model for the percolation is the random or 
ordinary percolation in which a vacant site or a vacant bond of the background lattice is 
randomly chosen to be occupied. The percolation transition in the random percolation is 
normally known to be continuous The percolation has been extensively studied during 
last 3 or 4 decades to be considered as a mature branch of sciences. 

However anomalous physical properties of exotic models on the percolation recently trig- 
gered some new studies. One kind of studies j^J was on the explosive percolation, which was 
first known to show supposedly discontinuous transition on the complete graph (CG) 2|, |3|. 
But subsequent studies on the explosive percolation have shown that the^ransition of the 
explosive percolation on CG or the mean-field transition is continuous 

Another kind of studies was on the agglomerative percolation (AP) 
cluster is randomly selected instead of a bond or a site. Then the selected cluster merges 
all the nearest neighboring clusters to form a new cluster. The phase transition in AP is 
shown to be continuous, but belongs to a new universality class different from the class 
of the random percolation if the base structure of AP is bipartite 121 ] . On the bipartite 
structure like a two-dimensional square lattice the merging process spontaneously breaks the 
Z 2 symmetry at the transition threshold, which is the origin of the new universality class 
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12J. In AP one 



12] . In contrast, the universality of the transition of AP on the triangular lattice, which is 



not bipartite, is the same as that of the random percolation l^ . Using analytical methods 
and numerical simulations, APs on the one-dimensional ring 8] , the two-dimensional square 
lattice and triangular lattice js], critical tree 10], and complex network 11 ] were studied. 
Through these studies AP on bipartite graphs is shown to belong to a new universality class 
different from that of the random percolation. 

To understand and establish a new universality class of the critical phenomena clearly and 
precisely the exact mean-field theory for the new model must be first understood. However 
the mean-field theory (MFT) for AP on bipartite networks was not clearly understood yet. 
To get MFT of AP, the analytic theory based on the generating function of the Erdos-Renyi 



(ER) random network was attempted 111 ] . However this analytic approach predicted the 
critical exponent 7 as 7 = 1/2, but from the numerical simulation on ER graph 7 = 0.88(10) 
is obtained, which is significantly larger than 7 = 1/2. This suggests a possibility that 



the analytic approach based on the generating function [ll| is still far from completion. 
Furthermore ER graph is not exactly bipartite. The numerical simulation study on the 
exact bipartite random graph earned only the critical exponent v and the fractal dimension 
D of the giant cluster as v = 4.7(2) and D = 0.567(6), which are close to those for ER 
network but differ by more than one standard deviation 12j. Therefore, at the present 
stage, MFT for AP on the bipartite graphs are far from completion. 

Recently the complete graph is widely used as a testbed for MFT {2], 4-6]. However the 
complete graph is not bipartite and one growth step of AP on the complete graph makes 
the entire graph a new single cluster. In contrast the Bethe lattice (infinite homogeneous 
Cayley tree) is the exact bipartite graph on which AP can be well defined. Moreover, 
the Bethe lattice is physically a very important substrate or medium on which MFTs for 
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various physical models become exact [14j . The ana 
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ic treatments of magnetic models |15| . 



16| on the Bethe lattice give important 



percolation (l|, |14(, localization [14j, and diffusion 
physical insights to subsequent developments of the corresponding research fields. Therefore, 
if the critical phenomena of AP on the Bethe lattice is completely understood, one knows 
MFT for AP exactly. 

One of the theoretical merits of the Bethe lattice is that one can set up exact self- 
consistent equations on the lattice. Recently we have developed an exact self-consistent 
simulation method for an arbitrary percolation process on the Bethe lattice Q]. From the 
self-consistent simulation method, we have shown that the phase transition of the Achlioptas- 
type explosive percolation [2| undergoes continuous transition regardless of the details of 
growth rules. In this paper we analyze the critical properties of AP on the Bethe lattice by 
use of the developed self-consistent simulation. In the self-consistent simulation the order 
parameter and the average size S of finite cluster on the Bethe lattice are directly 
measured. Therefore the exponents /3 and 7 are also obtained directly without the finite 
size scaling, and our work can indeed establish exact MFT of AP. 

In addition Lau et al. suggested the modified AP on the /c-partite graph, which we call 
APfc 12|. So AP2 means the original AP on the bipartite graph. Based on the simple 
arguments, the transition of AP fc is conjectured to belong to another new universality class 



when the growth of clusters in AP& spontaneously breaks the symmetry of the fc-partite 



graph 



12|. However the conjecture have never been confirmed quantitatively, yet. Therefore, 



in this paper we also study the MFT of AP 3 by using the triangular cactus structure, which 
is an expanded structure of the Bethe lattice and exactly tripartite [13I . By the self- 
consistent simulation we will also find the mean- field exponents /3 and 7 for AP 3 , or (3 3 and 
73. Finally from the results of AP2 from the Bethe lattice and AP3 from the triangular 
cactus, the inequalities between (/3 for AP2) and ^3 and between 72 (7 for AP2) and 73 
will be provided in the mean-field level. From the obtained inequalities, we will suggest the 
inequalities for all /Vs and 7^'s. 

This paper is organized as follows. The ordinary AP or AP2 on the Bethe lattice is 
studied based on the self-consistent simulation in Sec. II. AP3 on the triangular cactus is 
defined and studied Sec. III. Finally we summarize our results in Sec. IV. 



II. AP 2 ON THE BETHE LATTICE 




FIG. 1. Schematic diagram for AP2 on the the Bethe lattice with z = 3. The center part consists 
of a three- generation Cayley tree with edge sites denoted by ®. Each edge site is connected to an 
infinite cluster (IC) with the probability A or to a finite cluster (FC) of average size of Sf, with 
the probability 1 — A. Thick lines mean occupied bonds and thin lines mean vacant bonds. If the 
cluster B is selected, gray sites merge into one cluster. 

In the ordinary AP or AP 2 one cluster is randomly selected instead of a bond or a site and 
the cluster merges all the nearest neighboring clusters to form a new cluster. This means 
that in each growth process multiple bonds can be occupied at the same time. Therefore the 
natural control parameter in AP 2 is the number of clusters per site n instead of the fraction 
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p of the occupied bonds or sites 

The Bethe lattice is the infinite Cayley tree in which tree structures connected to the 
center site O are identical to one another as schematically shown in Fig. [1] The 
Bethe lattice is of course bipartite. Therefore AP2 on the Bethe lattice is expected to belong 
to a new universality class different from that of the random percolation. 

Let us now briefly explain the self-consistent simulation method for arbitrary percolation 
on the Bethe lattice with coordination number z in Ref. In the method we originally 
used the fraction p of occupied bonds or sites, but we use the number n of clusters per site 
in this paper. To set up self-consistent equations on which the simulation method is based, 
first consider a part of the Bethe lattice with m generations from the center site O, which 
has total N = 1 + z(k m — l)/(k — 1) sites, where k = z — 1. To make a complete Bethe 
lattice, one should add an infinite branch to each of zk m ~ l edge sites. To calculate the order 
parameter Poo{n) of percolation, which is defined by the probability for O to belong to an 
infinite cluster at a given n, we need to know the probability A with which an occupied edge 
site connected to an infinite cluster. Let P moo (n, A) be P^ which is calculated from a Bethe 
lattice with the m generations from O and zk m ~ x infinite branches. Then the self-consistent 
equation for P^ becomes 

-Poo = Pmoo{ n i A) = P m 'oo{ n i A) (1) 

for any combination of {m, m'}. Let us define P ms t{ n , A, S b ) as the probability that a cluster 
including O with s sites and t edge sites occurs within the m-generation tree. Then 

P moD (n, A) = 1 - ^(1 - Af Pmst(n, A, S b ) (2) 

t s 

where S b is the average size of the finite cluster connected to an edge site of the m-generation 
tree as shown in Fig. 1. The self-consistent equation for the average size S of the finite 
clusters including O can also be written as 

S = S m (n, A, S b ) = S m /(n, A, S b ) (3) 

where 

c , , <n E s ,tPmst(n,A,S b )[s + tS b }(l-AY 

S m (n, A, S b ) = ; — . (4) 

J- P OO 

If one cannot calculate P ms t(n, A, S b ) analytically to solve the self-consistent equations, one 
should estimate P ms t(n,A,S b ) indirectly. One of such indirect methods is a simulation 



method. We have developed a simulation method to solve self-consistent equations, which 
we call the self-consistent simulation In the self-consistent simulation, P mst {n, A, Sb) is 
estimated by the relation P mst (n, A, S b ) = N mst (n,A,S b )/N duster , where N mst (n,A,S b ) is 
the number of clusters including O with s sites and t edge sites within the m-generation 
tree that occurred in simulations. Of course, N c i uster is the total number of clusters which 
includes O within the m-generation tree that occurred in the same simulation runs. In the 
simulation both P mst (n, A, Sb) and P m > st (n, A, Sb) are estimated simultaneously using the 
Bethe lattice with the m-generation tree if m > ml . 

Since we don't know A and S b a priori, the iteration processes are needed in self-consistent 
simulation. From initially guessed values for A and S b , the final or saturated values of A 
and Sb are obtained by the iteration of unit simulation process. The unit simulation process 
consists of the following two steps. (I) By use of the simulation runs based on the given 
values A and S b , P mst (n, A, S b ) and P m 'st (n, A, S b ) are estimated. (II) From the estimated 
Pmst(n, A, S b ), new A and S b are calculated by utilizing self-consistent equations (CQ) and 
(J3J). In the unit simulation process to get new A and S b , the quantities like P mst (n, A, S b ) 
are estimated by averaging over at least 10 6 simulation runs. Such unit process is repeated 
until A and S b reach the saturation values. Using the saturated values of A and S b , P^ 
and S are estimated from Eqs. O) and (HI). In the self-consistent simulation, it should 

n 

be careful to choose m'(< m) for a given m as addressed in Ref. [7J. If m is too small, 
the clusters within the m'-generation tree cannot have physical properties of AP2 enough 
to give physically plausible solutions for self-consistent equations (JT]) and (j3J). If m! is very 
close to m, P moo (p,A) is numerically not so much distinct from P m > 00 {p,A) and the self- 
consistent equation (TjQ) hardly gives the physically right solution. From the simulations with 
various sets of {m, m'} it is confirmed that suitable choice of m! should be in the interval 
m/3 < m' < m/2. 

The results of the self-consistent simulation with z = 3, m = 20, m' — 9 for AP 2 are 
displayed in Figs. [2] and [31 From the data for P^ in Fig. [2] the critical density n c and the 
order parameter exponent (3 2 are obtained based on the equation 

Poo-K-n)^ 2 , (5) 

which holds for the ordered phase or for n < n c near the critical point, i.e., n — > n~ . The 
obtained n c and /3 2 are n c = 0.7223(1) and (3 2 = 1.79(3). We have checked the results for the 
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FIG. 2. (Color online) Plot of Poo for AP2 against n c — n with n c = 0.7223. The line denotes the 
relation (|5|) with n c = 0.7223 and 02 = 1.79. Inset is the raw-data plot of Poo against n. 

simulation for some other combinations (m = 14, m') with m/3 < m' < m/2 and found the 
same results. For another consistent checks we also applied the self-consistent simulation 
on the Bethe lattice with z = 6 to obtain 02 = 1.79(3). These numerical results for 02 is 
close to the previous estimate P2 = 1.78(8) on the ER graph [ll], but our estimate has much 
smaller errors. 




FIG. 3. (Color online) Plot of 5 for AP2 against \n c — n\ for ordered phase (n < n c ) and disordered 
phase (n > n c ). The lines denote the relations ([6]) with 7^ = 7^ = 0.88 and n c = 0.7223. Inset is 
the raw-data plot of S against n. 
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From the data for the average size S of finite clusters in Fig. [3] and the equation 

{In — n r |~ 7 2 if n r < n 
+ r ' (6) 

\n — n c \ 72 if n c > n 

we also estimated n c , 7^ and 73 ■ Obtained n c is nearly the same as that obtained from 
the data in Fig. [2j We also obtain 72 = 72" = 72" = 0.88(1), in which no asymmetry is 
found between the disordered phase (n c < n) and the ordered phase n c > n. The result 
72 = 0.88(1) is also consistent with the previous estimate 72 = 0.88(10) on the ER Graph. 
We also obtained 72 = 0.88(3) on the Bethe lattice with z — 6. These results clearly show 
that the obtained values of 02 and 72 for AP 2 are significantly different from those for the 
random percolation [l|. 

In conclusion our estimates 02 = 1-79(3) and 72 = 0.88(1) on the Bethe lattice are far 
more precise MFT exponents for AP 2 on the bipartite graph, since the dimensionality of the 
Bethe lattice is infinite. 



III. AP S ON THE TRIANGULAR CACTUS 



The triangular cactus was first introduced by Fisher and Essam 13] to investigate the 



effects of loops 14j on the percolation. As shown in Fig. H] (a), the triangular cactus with 
coordination number z = 4 can be constructed from the Bethe lattice with z = 3. Each site 
in the Bethe lattice is replaced with a triangle of three sites to form the triangular cactus 
as shown in Fig. [4] (a). Thus the dimensionality of the triangular cactus is infinite as the 
Bethe lattice. Moreover, the triangular cactus is exactly tripartite, not bipartite as shown in 
Fig. H] (b). Therefore one can expect that the critical phenomena of AP2 on the triangular 
cactus belong to the random percolation universality class. 

Recently, Lau et al. suggested the modified AP on the /c-partite graph, which we called 
APfc [12]. It is conjectured that the universality class of AP& depends on k 12] . However AP& 
for k > 3 has never been quantitatively studied, yet. In this section AP3 on the triangular 
cactus is studied to obtain MFT of AP3. 

In AP3 on the triangular cactus one cluster is randomly selected and the cluster merges 
some of the nearest neighboring clusters into a new cluster, instead of all neighboring clusters 
in AP 2 on bipartite graph. In a tripartite graph, initially, three colors are arranged such that 
no pair of nearest neighbor sites has the same color. Therefore we can identify the cluster 
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FIG. 4. (Color online) (a) Formation of the triangular cactus from the Bethe lattice with z = 3. In 
the cactus each site of the Bethe lattice is replaced with a triangle composed of three sites. Each 
edge site which is denoted by "X" is connected to an infinite branch, (b) AP3 on the triangular 
cactus. The triangular cactus is tripartite as shown in the figure. If a blue-colored cluster indicated 
by an arrow is selected as in the left figure, it agglomerates all "R" neighbors and becomes a "R"- 
colored cluster by the rule, R — ^ G — > B — >■ R as in the right figure. 



by colors such as red (R), green (G), and blue (B). AP 3 is denned such that a selected 
cluster with "R" can join only with neighbors of the color "G", "G" can join only with 
neighbors of the color "B", "B" can join only with neighbors "R", based on a cyclic rule, 
R — > G — > B — > R 12|. For example, a certain cluster with "R" is selected, then the cluster 
merges all neighboring clusters colored by "G" into a new cluster, and the merged cluster 
becomes a new "G"-colored cluster from the rule, R — > G — > B — > R. One can apply the 
R — y B — 7- G — > R rule to the model, but it cannot be physically different from the model 
with the R — > G — > B — > R rule. 

For MFT of AP 3 on the tripartite graph, we use the self-consistent simulation method for 
AP3 on the triangular cactus. The self-consistent simulation is almost the same as that for 
AP2 on the Bethe lattice. First consider a part of the triangular cactus with m-generations 
from O, which has total number of sites N = 2 m+2 — 3. To make a complete triangular 
cactus, one should add an infinite branch to each of 2 m+1 edge sites. Other details of the 

9 



self-consistent simulation on the triangular cactus are exactly the same as those on the Bethe 
lattice. For instance, Eqs. (JT|)— (SJ) are the self-consistent equations not only for the Bethe 
lattice but for the triangular cactus. 

We first checked the results of self-consistent simulation with m = 20, m' — 9 for AP 
on the triangular cactus. From the data for P^ and S, n c , (3 and 7 is estimated as n c = 
0.6761(1), (3 = 1.01(2), and 7 = 1.00(2). Since the triangular cactus is not bipartite but 
tripartite, this result supports that the critical phenomena of ordinary AP on the tripartite 
graph belong to the random percolation universality class with (3 = 7 = 1 as expected in 



Ref. 



12|. 
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FIG. 5. (Color on line) Plot of for AP3 on the triangular cactus against n c — n with n c = 0.5681. 
The line denotes the relation similar to ([5]) with n c = 0.7223 and /?3 = 1.35. Inset is the raw-data 
plot of Pqo against n. 

On the other hand, the critical phenomena of AP3 is different from AP on the triangular 
cactus. The results of the same self-consistent simulation for AP 3 on the triangular cactus 
are shown in Figs. |5] and [61 By using the similar equations to Eqs. ()5]) and ([6]) we have 
obtained the order parameter exponent $3 and the susceptibility exponent 73 for AP3 on 
the triangular cactus. The results are (3 3 = 1.35(5) and 73 = 0.94(2) with n c = 0.5681(1) as 
shown in Figs. [5] and [6J In conclusion our estimated [3% and 73 are the first MFT exponents 
for AP3, since the dimensionality of the triangular cactus is infinite. 

The obtained exponents (3 3 and 73 satisfy the relations (3^ < (3 3 < (3 2 and 7^ > 73 > 72, 
where /3<x>( = 1) an d 7oo(= 1) are the MFT exponents of the random percolation. The 
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FIG. 6. (Color on line) Plot of S for AP3 on the triangular cactus against n c — n with with 
n c = 0.5681 for disordered phase (n > n c ). S for ordered phase (n < n c ) is not shown, because it 
behaves almost the same as that for disordered phase. The line denotes the relation similar to Eq. 
([6]) with 7^ = 0.94(2). Inset is the raw-data plot of S against n. 

relations of exponents suggests that the MFT exponents of AP& on the fc-partite graph 
approach to those of the random percolation as k increases. 



IV. SUMMARY 

Finding the exact MFT is the first step to understand the various physical properties 
of a new model. AP was suggested as a natural extension of the standard percolation 
model. Some numerical studied for AP have been done on lower-dimensional lattices and 
random graphs [si, fiol 12] . Based on those numerical studies, it was conjectured that AP 



belongs to a new universality class if the growth of cluster breaks symmetry on /c-partited 
graph. However, the mean-field approach based on the evolutionary dynamics of clusters 
did not agree with the numerical simulations. This strongly indicates that AP is not fully 
understood even in mean- field level 111. Therefore, in order to provide an exact MFT, we 

Q 

apply the self-consistent simulation method [7fl to APs on the Bethe lattice and triangular 
cactus. From the direct and precise measurement of and 5* through the self-consistent 
simulation, we obtain (3 2 = 1.79(3) and 7 2 = 0.88(1) on the Bethe lattice when Z 2 symmetry 
is broken spontaneously at the transition threshold. Similarly, we obtain /3 3 = 1.35(5) and 
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73 = 0.94(2) on triangular cactus if Z3 symmetry is broken spontaneously. However, since 
the triangular cactus is not bipartite, AP model on triangular cactus gives (3 = 1.01(2) and 
7 = 1.00(2). This result shows that ordinary AP on triangular cactus belongs to the same 
universality class with ordinary percolation. Therefore, the results for AP 3 on triangular 
cactus provide the exact MFT verifying the Lau et a/.'s arguments [3]. In addition, by 
comparing the obtained critical exponents with those of random percolation, we also find 
the inequalities (3^ < (3 3 < (3 2 and 700 > 73 > 72- These inequalities also quantitatively 
verify the conjecture that the universality class of AP& depends on k 121 ] . 
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